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A LAGRANGIAN APPROACH FOR THE INCOMPRESSIBLE 
NAVIER-STOKES EQUATIONS WITH VARIABLE DENSITY 

RAPHAEL DANCHIN AND PIOTR BOGUSLAW MUCHA 

Abstract. Here we investigate the Cauchy problem for the inhomogeneous Navier-Stokes equa- 
tions in the whole n-dimensional space. Under some smallness assumption on the data, we show 
the existence of global-in-time unique solutions in a critical functional framework. The initial 
density is required to belong to the multiplier space of -B"/ ~ (R"). In particular, piecewise 
constant initial densities are admissible data provided the jump at the interface is small enough, 
and generate global unique solutions with piecewise constant densities. Using Lagrangian co- 
ordinates is the key to our results as it enables us to solve the system by means of the basic 
contraction mapping theorem. As a consequence, conditions for uniqueness are the same as for 
existence. 
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Introduction 
We address the well-posedness issue for the the incompressible Navier-Stokes equations with 



> 

\^ I variable density in the whole space 



O ■ (0-1) 



^ 



dtp + u ■ V p = 0, 

p{dtu + u ■ Vu) - fiAu + VP = 0, 

divn = 0, 

u\\t=Q = Uo. 

Above p = p{t,x) € M+ stands for the density, u = u{t,x) € M^, for the velocity field and 
P = P{t, x) G M, for the pressure field. The viscosity coefficient /i is a given positive real 
cd '_ number. We supplement this system with the following boundary conditions: 

• the velocity u tends to at infinity, 

• the density tends to some positive constant p* at infinity. 

The exact meaning of those boundary conditions will be given by the functional framework in 
which we shall solve the system. In what follows, we will take p* = 1 to simplify the presentation. 

This old system has known a renewed interest recently, among the mathematics community. 
The existence of strong smooth solutions with positive density has been established in e.g. [11] 
whereas the theory of global weak solutions with finite energy has been performed in the book 
[12] by P.-L. Lions (see also the references therein, and the monograph [3]). As pointed out in 
[5], it is possible to construct strong unique solutions for some classes of smooth enough data 
with vanishing density. 

In the present paper, we aim at solving the above system in critical functional spaces, that is 
in spaces which have the same invariance with respect to time and space dilation as the system 
itself (see e.g. [6] or [7] for more explanations about this nowadays classical approach). In this 
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2 R. DANCHIN AND P.B. MUCHA 

framework, it has been stated in [1, 6] that, for data {pQ,uo) such that 

(po - 1) € ^;//(M"), uo G ^^.(''"^(K") with divno = 

and that, for a small enough constant c, 

(0-2) IIPO - l|l^"/P(IR>n) + ^""^ll^oll^n/p-1 < C, 

we have for any p G [1, 2n) : 

• existence of a global solution (/?, n, VP) with p — 1 € Chm+;Bp-^), u^Chi^+^B^i ) 

and dtu, V^u, VP G Li(M+; P^/i^"^) ; 

• uniqueness in the above space if in addition p < n. 

Those results have been somewhat extended in [2] where it has been noticed that po — 1 may 
be taken in a larger Besov space, with another Lebesgue exponent. 

The above results are based on maximal regularity estimates in Besov spaces for the evolu- 
tionary Stokes system, and on the Schauder-Tychonoff fixed point theorem. In effect, owing 
to the hyperbolicity of the density equation, there is a loss of one derivative in the stability 
estimates thus precluding the use of the contraction mapping (or Banach fixed point) theorem. 
As a consequence, the conditions for uniqueness are stronger than those for existence. 

Let us also point out that all the above results concerning existence with uniqueness require 
the density to be at least uniformly continuous. This condition has been somewhat weakened 
recently by P. Germain in [9]. However, there, initial densities with jump across an interface 
cannot be considered. 

In the present paper, we aim at solving System (0.1) in the Lagrangian coordinates. The 
main motivation is that the density is constant along the flow so that only the (parabolic type) 
equation for the velocity has to be considered. We shall show that, after performing this change 
of coordinates, solving (0.1) may be done by means of the Banach fixed point theorem. As a 
consequence, the condition for uniqueness need not be stronger than that for the existence, and 
the flow map is Lipschitz continuous. 

Our main result states the global in time existence of regular solutions to the inhomogeneous 
Navier-Stokes equations in M" in optimal Besov setting, under suitable smallness of the data. As 
regards the initial density, the admissible regularity is so low it may have (small) jumps across 
a C^ interface. This is of particular interest from the viewpoint of physics as it implies that 
motion of a mixture of two incompressible fluids with slightly different densities can be modeled 
by the inhomogeneous Navier-Stokes equations. In addition, the regularity of the constructed 
velocity suffices to preserve the C^ regularity of the interface between fluids. 

We now come to the plan of the paper. In the next section, we present our main results and 
give some insight of the proof. Section 2 is devoted to solving the linearized system (0.1) in 
Lagrangian coordinates. This will enable us to define a map ^ : Ep ^>- Ep where Ep stands 
for the functional space in which the Lagrangian version of the momentum equation of (0.1) 
is going to be solved. That <I> fulfills the conditions of the contraction mapping theorem on a 
small enough ball of Ep is the main purpose of Section 3. In the Appendix we prove several 
important results concerning the Lagrangian coordinates and Besov spaces. 

Notation: Throughout, the notation C stands for a generic constant (the meaning of which 
depends on the context), and we sometimes write X <Y instead oi X < CY. Finally, for 
A = {Aij)i<ij<n and B = {Bij)i<ij<n two matrices, we denote A : B = ^-j AijBji. 
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1. Main results and principle of the proof 

Let us first derive formally the Lagrangian equations corresponding to (0.1) . Let X„ be the 
flow associated to the vector-field u, that is the solution to 

(LI) Xuit, y) = v+ I u{t, Xu{t, y)) dr. 

Jo 

Let us recall that by Liouville's formula for transport equations, the divergence-free condition 

is equivalent to |-DXu| = 1. In other words the map (1.1) is measure preserving. Now, denoting 

p{t,y):=p{t,Xu{t,y)), P{t,y) := P{t,Xu{t,y)) and u{t,y) = u{t,Xu{t,y)) 

with (pjUj'VP) a solution of (0.1), and using the chain rule and Lemma 1 from the Appendix, 
we gather that p{t, •) = po and that (w, V-P) satisfies 

j podtu-pdivy{A,jAuVyu) + ^AuVyP = with Au = {DyXuy^, 
[ divj^(^„n) = 0. 

Motivated by prior works (see e.g. [1, 2, 6, 7]) we want to solve the above system in critical 
homogeneous Besov spaces. Let us recall that, for 1 < p < oo and s < n/p, a tempered 
distribution u over M" belongs to the homogeneous Besov space i?p]^(]R") if 

u = Y,^jU in 5'(M") 
and 

Here {Aj)j^z denotes a homogeneous dyadic resolution of unity in Fourier variables (see e.g. 
[4], Chap. 2 for more details). 

Loosely speaking, a function belongs to Bpi{M."') if it as s derivatives in Lp(]R"). In the 
present paper, we shall make an extensive use of the following classical properties: 

• the Besov space n'^i{M.'^) is a Banach algebra embedded in the set of continuous func- 
tions going to at infinity, whenever 1 < p < oo; 

• the usual product maps 5p(*'"\]R") x ^^//(M") in ^^(''"^R") whenever 1 < p < 2n. 

From now on, we shall omit M" in the notation for Besov spaces. We shall obtain the existence 
and uniqueness of a global solution (u, VP) for (1.2) in the space 

Ep := {{u,VP) : u € a(M+; ^;//-'), dtU,\/\VP € Li(M+; ^;/^-')}, 
and we shall endow Ep with the norm 

\\{u,VP)\\e, := ML^(^^^,B;/r') + ll^*^'^^'^'^^llLi(R+;B;/r')- 

We shall also use the local version Ep{T) of Ep, pertaining to functions defined on [0,T) x M". 
Writing out the exact definition and the corresponding norm is left to the reader. 

The required regularity for the initial density pQ is that it belongs to the multiplier space 
M.[Bp^ ) for B^i , that is the set of those distributions pQ such that il'Po is in B^^ 
whenever ip is in -B" / , endowed with the norm 
(1-3) \\p^\\m{b"'^-^) •= supllV'Poll^n/p-i 

where the supremum is taken over those functions ip in B^i with norm 1. 



The reader may refer to the Appendix for the rigorous derivation in our functional setting. 
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Let US now state our main result. 

Theorem 1. Let p € [l,2n) and uq be a divergence- free vector field in B^i (M"). Assume 

that the initial density po belongs to the multiplier space M.[B^'^ ). There exists a constant c 
depending only on p and on n such that if 

then System (1.2) has a unique global solution (u, VP) in Ep. Moreover, we have 



iu,VP)<C\\uo\ 



B, 



n/p— 1 

p,i 



for some constant C depending only on n and on p, and the flow map {po,uq) i — > (n, VP) is 
Lipschitz continuous from A4{B^i ) x -B"/ to Ep. 

In the case where only the density satisfies the smallness condition, we get the following 
local-in-time existence result: 

Theorem 2. Under the above regularity assumptions, there exists a constant c depending only 
on p and on n such that if 

then there exists some T > such that System (1.2) has a unique local solution (u, VP) in 
Ep{T), and the flow map [pQ, uq) i — > [u, VP) is Lipschitz continuous from M.{B'^J^ ) x -B" / 
to Ep{T). 

The regularity given by Theorem 1 ensures that the map defined in (1.1) is defined globally 
(see the Appendix). Coming back to the Eulerian formulation, this will enable us to get the 
following result : 

Theorem 3. Under the above assumptions. System (0.1) has a unique global solution (/?, n, VP) 
with p € Loo(M+;M(Pp//"^)) and (n, VP) G Ep. 

Let us make a few comments concerning the above assumptions. 

• The condition 1 < p < 2n is a consequence of the product laws in Besov spaces. Let us 
emphasize that any space L^ n Bg^r^ with q satisfying 

. X 111 ,111 

1.4 - > and - > , 

q n p q p n 

embeds in A4{B^f ), a consequence of basic continuity results for the paraproduct 
operator (see [4]). Hence the above statement improves those of [1, 2] as regards the 
uniqueness. In particular, one may take the initial velocity in a Besov space with a 
negative index of regularity, so that a highly oscillating "large" velocity may give rise to 
a unique global solution. 

• In contrast with the results of [1, 2, 6] it is not clear that the above statements may be 
generalized as so if the viscosity depends on the density. Recall that in this case, the 
diffusion term in the momentum equation of (0.1) reads div {p{p)(\7u + Vn)) where 
/U is a given suitably smooth nonnegative function. One may easily extend the above 
statements to this case under the stronger condition that {p{po) — m(1)) is small in 

M(p;//). 



We here consider only the case of small data to simplify the presentation. 
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The space A4{Bpf ) is in fact much larger than L^o ^ B^'oi with q satisfying (1.4) (see 
e.g. [13], Chap. 4). It contains characteristic functions of C^ bounded domains, whenever 
p > n — 1 (see the proof in Lemma 5). Hence, our result applies to mixture of fluids, which 
is of course of great physical interest. In addition, given that the constructed velocity field u 
is divergence- free and admits a C^ flow X (see again the appendix), we deduce the following 
result which emphasizes the range of Theorem 1 (we just state the case of small velocities to 
simplify the presentation): 

Corollary 1. Assume that uq € -Bpi with divuo = and n — 1 < p < 2n. Let ilo be a 
hounded C^ domain of M". There exist two constants c (depending only on p and n) and d 
(depending only on p, n and VLq) such that if 

(1.5) IkolLn/p-i < c and po = 1 + c^Xno with \a\ < c', 

then System (0.1) has a unique global solution {p^u^VP) with p G Lqo (IK+ ; A4 (i?" j^^ )) and 
{u, VP) S Ep. In addition, for all time, 

(1.6) p{t) = 1 + axut where Qt = Xu{t, ^o). 
Besides, the measure and the C^ regularity of di^t o-t^^ preserved for all time. 

Let us give the main ideas of the proof of existence. Obviously, it suffices to find a fixed point 
for the map : [q, v) i— > (/?, u) where (/?, u) stands for the solution to the linear system 

dtp + v ■Vp = 0, 

q{dtu + V ■ Vn) - pAu + VP = 0, 

divn = 0, 



(1.7) 



U\lt=o = Uq. 

Although it is possible to prove uniform estimates in Loo(K+; A^(Pp i )) x Ep for {p,u,VP), 
we do not know how to get stability estimates in the same space, owing to the hyperbolic nature 
of the density equation. As a consequence, the contraction mapping theorem does not apply. 

In the present paper, we shall rather define the solution of the above system in the Lagrangian 
coordinates corresponding to v. For such coordinates, the density is time-independent. So, given 
some reference vector field v and pressure field VQ with {v, VQ) G Ep^ one may define (u, VP) 
to be the solution in the Lagrangian coordinates y = X^'^{t,x) pertaining to v of the linear 
system (1.7) (that X^ is a C"*^ -diffeomorphism over M" is proved in the appendix). 

Let us give more details. We assume that |-DX^| = 1 and set 

p{t,y):=p{t,X^{t,y)), P{t,y) := P{t,X.,{t,y)) and u{t,y) = u{t,X^{t,y)), 

where {p,u,VP) stands for a solution to (1.7). Then we have p{t, •) = po ^-^d (see the proof in 
appendix) 

Podtu - pdw {AjA^Vyu) + ^A^VyP = 0, 
(1.8) 

divj^(y4^u) = 

with 

(1.9) A^ = {DyX^y^ and X^{t,y) = y + [ v{T,y) dr. 

Jo 

Solving this linear system globally turns out to be possible under some smallness condition 

over V and po — 1. This will enable us to define a self-map $ : (vj'VQ) i— )■ (u, VP) on Ep. 

Then it will be only a matter of checking that if the data po ^-nd ^o satisfy a suitable smallness 

condition then the map <I> fulfills the assumptions of the standard Banach fixed point theorem. 
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The key to that will be estimates for the Stokes system in Ep (see Proposition 1) and a "magic" 
algebraic relation involving the second equation of (1.8) (which has been used before in e.g. 
[14, 15] in a different context). 

2. The linear theory 

Our proof of existence for the linear System (1.8) will be based on the following a priori 
estimates for the Stokes system, the proof of which may be found in [7] : 

Proposition 1. Let p € [l,oo] and s £ R. Let uq € ^^,i(M") and f £ Li(0,T; ^^ ^(M")). Let 
5 : [0, T] X M" ^ M be such that 

V5eLi(0,r;^^i(IR")), dtg = d[vR with R e Li{0,T; B^-^iR'')) 
and that the compatibility condition g\t=o = divuQ on M" is satisfied. 
Then System 

dtu - fiAu + VP = f in (0, T) x M" 

(2.1) divu = g in (0,T)xIR" 

u\t=to = uo on R" 

has a unique solution {u, VP) with 

u£C{[0,T);B'p^^{W)) and dtu,V^u,VP e Li{0,T; B^^^{W)) 
and the following estimate is valid: 

(2.2) < C'dl/, ^V^, i?||^^(-Q^j,.^s^(I[j„)) + ll'"0|lBs^(]Rn)) 

where C is an absolute constant with no dependence on fi and T. 

Granted with the above statement, we rewrite System (1.8) as 

dtu - fiAu + VP = il-po)dtu + ^div {{AjA^ - Id)Vn) + (Id - ^A^)VP 

(2.3) div u = div ((Id - A^)u) 
u\t=o = Uq. 

We assume that the vector-field v from which A^ and DX^ are defined satisfies 

(2.4) veC{R+;B;fr^), dtv,V^v G Li{R+;B;lr^), \DX,\ = 1 
and that, for a small enough constant c. 



(2.5) / \\Dv\\-„/pdt<c. 

Jo P'l 

Even though this system is linear, it cannot be solved directly by means of Proposition 1 for 
the right-hand side depends on the solution itself. So in order to prove the existence of -u, we 
shall look for a fixed point of the map 

^:(tD,VQ)^(n,VP) 



Because homogeneous 0-th order muhipUers act on any homogeneous Besov spaces Bp i, one may take any 
indices s and exponents p £ [l,oo], in Proposition 1. 
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where (u), VQ) € Ep and (n, VP) stands for the solution of 

' dtu - fiAu + VP = fiw, VQ), 
(2-6) < div-u = g{w), 

, u\t=o = Uq. 
Above, gi^w) := div ((Id — Ay)w) and 

(2.7) fiw, VQ) := (1 - po)dtw + ^div {{AJA - M)Vw) + (Id - ^A^)VQ. 

We claim that, if v satisfies (2.4) and the smallness condition (2.5), then for any {w,VQ) in 
the space Ep with 1 < p < 2n, the above system has a unique solution (u, VP) in Ep and that, 
in addition, the map ^ fulfills the required conditions for applying the contracting mapping 
theorem. 

The existence of (n, VP) will stem from Proposition 1 once it has been checked that f{w, VQ) 
and g{w) fulfill the required conditions. As regards g{w), this stems from the following "magic 
formula" 

(2.8) giw) = div ((Id - A^)w) = Did : (Id - A^), 

a consequence of the fact that |I?X^| = 1 (see Corollary 2 in the Appendix). 

Bounds for g{w) . Let us first check that g{w) € Li(M+; B'l^). As Div G Li(]R+; dIi) and as, 

according to (A. 8), Id — A^ is in Loo(IK+; d^i), this is a consequence of the fact that P" / is 
a Banach algebra and that 

g{iB) = Div:{M-A^). 

In addition, we get 

(2-9) \\g{iv)\\,,^^,^^/,<\\Dv\\^,^^.^/,.\\DiB 



n/ps 



Next, we see that dt{g{iB)) = divP^(it)) + divP^(ttj) with 

R^iiu) := (Id - A^)dtiu and R^{id) := -dtA^iu. 

So, according to (A. 7), (A. 9) and because the product operator maps B^-^ x P"/ in P"/ 
whenever p < 2n, we see that R^{iJu) and R^{iv) belong to Li{M.^; B^ -^ ) and that 

(2.10) ll^'(^)llLi(R+;B;_/r^) ^ ll^^llLi(R+;B;/-)ll^*^ilLi(R+;ij;//-i)' 

(2.11) ll^'(^)llL,(R+;B;,/r^) ^ "^^llLi(R+;B;.//)ll^llL^{R+;B;,/r^)- 

Bounds for f{iu, VQ) ■ That the first term of f{vj) belongs to Li(IR+; P" / ) is a consequence 
of the definition of the multiplier space M (Pp i ) ; in addition we have 

(2.12) \\{l- PQ)dtW\\^r./r,-l < 111 -poll... Rn/p-lJ|5tU) 



S^y-. ^ ,1^ ^"H^(B^y-) 



/p-1. 



n/p—1-^ 



Next, according to (A. 11), the second term of f{w) belongs to Li(M+;Pp^ ) and 

(2.13) lldiv {AjA, - Id)Vw))||^^(jj^.^„/^i^ < ll^^^lli^(K+.B;_/P)P^lli^(u+;B;/^)- 

Finally, Inequality (A. 8) and the fact that the product operator maps B^i x B^i in 
Bpf ii p < 2n ensure that 

(2.14) ||(Id - ^^.)VQ||,^(^^^^n/.-i) < l|V^II^^(«^^^./.)l|VQ||,^(^^^^./.-i). 
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So putting together (2.9) to (2.14), one may conclude from Proposition 1 that for any (w, VQ) 
in Ep^ System (2.6) has a unique solution [u,'S/P) in Ep. In addition 



\\{u,VP)\\e^ < C[\\Uq\\ ^r./p-l + (111 -P0ll_^(^n/P-1) + ll^^lli,(IR+;B"/P))ll(^'V<3)ll£;p 

As a consequence, there exists a positive constant c (depending only on n and on p) such 
that if (2.5) is satisfied and 

(2-15) I|l-P0||_;^(^n/P-1) <c, 

then 

(2.16) \\^{w,VQ)\\e, < C\\u4^^j,-, + ]^\\{w,VQ)\\e,. 

Banach theorem thus entails that the linear map ^ admits a unique fixed point in Ep^ that we 
shall still denote by {u^VP). Let us emphasize that Inequality (2.16) ensures that 

(2.17) ||(n,VP)|b, <2C||no||^„/p-i, 

p,i 

and that, by construction, div A^u = hence, according to Corollary 2, |DX„| = 1. In other 
words, given v fulfilling (2.4) and (2.5), System (2.3) admits a unique solution {u^VP) fulfilling 
the same conditions, together with (2.17). 

3. The inhomogeneous Navier-Stokes equations 
Let us denote by Ep the closed subset of Ep containing all the couples {v, VQ) such that 

\DX^\ = l and \\{v,VQ)\\e,<R. 

According to the previous section, if one takes (€», VQ) in Ep with v satisfying |-DX^| = 1 
and (2.5), then (1.8) admits a solution [u^VP) in the same space, such that |-DX„| = 1. Let 
^{v, VQ) denote this solution*^. We claim that if uq is small enough with respect to ji in i?" ^ > 
if the density po satisfies (2.15) and if R is small enough, then $ admits a unique fixed point 
in E^, as a consequence of the contracting mapping theorem. 

3.1. Stability of a small ball of Ep by <I>. Assume that (2.15) is satisfied and let us take 
R = cfi. Then (2.5) holds true whenever (vj'VQ) is in Ep. Therefore {u,VP) := ^{vjVQ) 
satisfies (2.17). So it is clear that if 

(3.1) 2C\\uo\\.n/p-i < cfi, 

with C as in (2.17), then (u, VP) is in E^, too. 

3.2. Contraction properties. In this part, we show that under Conditions (2.15) and (3.1) 
(with a greater constant C and smaller constant c if needed), the map <I> : E^ — )■ Ep is 
1/2-Lipschitz. 

So we are given (wi, VQi) and {v2,'VQ2) in Ep, and denote 

(ni, VPi) := <5(^i, VQi) and {u2, VP2) ■= Hv2, VQ2). 

Let Xi and X2 be the flows associated to vi and V2. Set Ai = {DXi)^"^ for i = 1,2. The 
equations satisfled by 5a := U2 — ui and V6P := VP2 — VPi read 

J dtSu - fiAdu + V6P = 6f := 6fi + 6/2 + Sh + /udiv 6fi + /idiv 5/5 

1 divdu = dg := div ((Id - A2)Su + {Ai - A2)ui) 



Of course, it is independent of VQ. However, prescribing the pressure is needed so as to define a map from 
a subset of Ep to itself. 
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with 

6fi := {1 - po)dtdu, (5/2 := (Id - %)Vff', <5/3 := Vi " ^2) VPi, 

Sf4 := {A2^A2 - AjA^) Vui, 5h ■= [A^^A^ - Id)V5u. 

Once again, bounding ((5u, VdP) will stem from Proposition 1, which ensures that, if dtbg = 
divdR then 

(3.2) \\{du,V5P)\\E, < 115/11^^ (j^^.^„/P-i) +^II<^I!^^(]R^.^;/P) + \\6R\\^^^^^,^„j,-iy 

So we have to estimate 6/1,6/2,6/3 and 6/4,6/5 in Li{M.-^-; B^-^ ) and Li{E.^;Bp^), respec- 
tively. First, from the definition of the multiplier space M.{B^J^ ), we readily have 

(3.3) ll'^/lllLi(R+;i?;//-^) ^ 11^0 - ^II.M(B;/r')"^*'^"ii (!«+;</"')• 

Next, using Inequalities (A. 8), (A. 11), and product laws in Besov space yields 

Inequality (A. 13) ensures that 

whereas Inequalities (A. 13), (A. 14) yield 

(^•^) \\^f^\\L,(R+;B;/^r) ^ II^*IIli(R+;B;/-)II^^iIIli(R+;</)- 

In order to bound 6g in Li(R+; B^-^), we shall use the fact that, by construction, 

divtij = div ((Id — ylj)uj) = Dui : (Id — A,i). 

Hence 

6g = D6u:{ld- A2) - Dui : {A2 - Ai). 
Now, easy computations based on (A. 8) and (A. 13) yield 

(3.8) \\D6u : (Id - ^2)||^^(j,^.^n/P) < \\Dv2\\^^^^^.^u/,^\\D6u\\^^^^^,^„/,y 

(3.9) \\Dui : iA2 - Ai)\\^^^^^,^^/,^ < \\Dui\\^^^^^.^^/,^\\D6v\\^^^^^,^^/,y 

Finally, to bound dt6g, we decompose it into div {6R1 + 6R2 + 6R3 + 6R4) with 
SR^ = -dtA2 6u, 6R2 = (Id - A2)dt6u, 

6R3 = dtiAi - A2) ui, 6R4 = (Ai - A2)dtui. 
Using (A. 7), (A. 9) and product laws in Besov spaces, we see that 

(^•^°) II'^-^iIIli(r+;</-^) ^ \\^^^k,iR+-,B;^r)^\^k^(R+;B;/r'y 

In order to bound 6Rs, it suffices to take advantage of (A. 15). We get 
Finally, using again (A. 14), we see that 

(3.13) \\^^^k,iR+;B;^r') ^ II^*IIli(M+;B;//) 11^*^1 llLi(R+;B;/r')- 
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One can now plug Inequalities (3.3) to (3.13) in (3.2). We end up with 

\\{du,V6P)\\E, < C(||l -po||_^(5n/.-i) + \\Dv2\\^^^^^.^^.^/r,^)\\{du,V6P)\\E, 

So we see that if (2.5) and (2.15) are satisfied for i'i,U2 and po with a small enough constant 
c, then we have 

\\{du,V6P)\\E, < 2CRf,-'\\{6v,VQ)\\E,). 

Hence, the map <^ : Ep i-^ Ep is 1/2-Lipschitz whenever R and the data have been chosen so 
that (2.15), (3.1) are satisfied and iCR < fi. This completes the proof of existence of a unique 
solution to System (1.2) in E^. 

3.3. Stability estimates in Ep. In this part, we want to prove stability estimates in Ep for 
the solutions to (1.2). This will ensure both uniqueness and that the flow map is Lipschitz. 

So we consider two initial divergence-free velocity fields uo,i and uo,2 in ^p i > and den- 
sities {po,i,Po,2) in ■^{Bpi ) satisfying (2.15) and (3.1). We want to compare two solutions 
(tti,V-Pi) and {u2.,VP2) in Ep of System (1.2), corresponding to data (/Oo,i)^o,i) and (/Oo,2 5 1*0,2) • 

The proof is similar to that of the contractivity of <5 : we have to bound 

aj{t) := !l'5u|li^(o,t;B;_/p-i) + \\dt5a,pV^5a,VmL,i^,,t;B;/r^) 
(with Su := U2 — ui and V5P := VP2 — VPi) in terms of ||&io|| o"/p-i and ||<^Po|| yy^/rj^/p-i-j- 
Now, the system for {5u, VdP) reads 

r dtdu - pASu + V5P = 5fo + 5fi + 5/2 + dfs + pdiv {5U + 6/5), 

\ div (5u, = div ((Id - A2)&i + (^1-^2)^1) = DSu : (Id - A2) + Dui : {A1-A2) 

with (5/o := 5pQ dtui and where J/j for i G {1, • • • , 5} has been defined in the previous subsection. 
Of course, now the matrices Ai and A2 correspond to the vector-fields ui and U2- 

Using Proposition 1, we gather that for all t G [0,T) 

3 5 4 

4=0 ' j=4 ' i=l 

(3.14) +fi\\Ddu : (Id-^2)||^^(o^,.^n/P) +^ll^i^i : (^2-^i)||^^(o_,.^n/P). 

From the definition of the multiplier space A4{B'^f ), we readily have 

lll^/olU"/?-! < W^PoWj^/ryri/p-lAldtUlW^n/p-l. 

The other terms may be bounded as in the previous subsection. So we eventually conclude that 
for aU t£[0,T) 

whenever, for i = 1,2, 

sup \\Ui{t)\\^n/p-l + (p^^\\dtUi,VPi\\^n/p-l + \\DUi\\^n/p]dt 

tG [o,r] p.i Jo ^ p.i p.i ^ 

is small enough. 

This latter condition is a consequence of (3.1). This completes the proof of stability estimates. 
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3.4. Proof of the local-in-time existence result. We here explain how the arguments of 
the previous subsections have to be modified so as to handle large initial velocities. 

Let us first notice that the computations that have been performed in Section 2 also hold 
locally on [0,T) whenever v satisfies 

(3.15) V € a([0,r);^;/*'-'), dtv,V^v £ Li(0,T; ^;//"'), \DX,\ = 1 on [0,T) x M" 

and 

rT 



(3.16) / 

Jo 



Dv\\^„/p dt < c. 
p.i 



This ensures that, under Condition (2.15), System (1.8) may be solved locally in Ep{T). Of 
course. Inequality (2.17) is still satisfied (for the norm in Ep(T)). However, it is not accurate 
enough so as to solve the nonlinear system, if uq is too large. To overcome this, we shall apply 
the contracting mapping theorem in some suitable neighborhood of the solution {ul,VPl) to 
the "free" Stokes system, that is 

ajUL - ^AuL + VPl = in [0,r)xIR", 

(3.17) divnL = in [0,r)xIR", 

UL\t=o = uo on M". 

Setting (n, VP) := ^{v^VQ) , we want to show that if T is small enough (a condition which 
will be expressed in terms of the free solution only) then {u,VP) := (n — ul,V{P — Pl)) is 
small. For that, we shall apply Proposition 1 to the system satisfied by {u,VP), namely 

dtu - fiAu + VP = /(n, VP) 

(3.18) divu = g{u), 

u\t=o = 

where f{u,VP) and g{u) have been defined in Section 2. 

On the one hand, we shall bound f{u,VP) in Li{0,T; Bpf ) and g{u) in Li{0,T; Bpf) 
exactly as in Section 2, on the other hand, decomposing dt{g{u)) into R^{u) + P?{u), we see 
that the bound (2.11) for E?'{u) is not accurate enough as it involves ||n|L ,„t, 6"/p-l^ which 

need not be small for T going to 0, if uq is large. So we shall rather write 

R^{u) = -dtA^ UL - dtA^ u, 
and use product laws and Inequality (A. 10), to get 

So, finally, using Proposition 1 and decomposing everywhere u and VP in ul + u and 
VPl + VP, we get 



l,DP)\\E^iT) < C(||l -/>0||_^(^n/p-l) + ||i^l^||^^(o^y.^n/P)) 



x{\\{u,DP)\\e^(t) + \\dtUL, fiD'^UL, DPL\\^^^Qr^,j^r./p-i^) + C'||i'||^^(Q^y.^n/p)||nL||^^^g^^.^„/p^. 

Therefore, if (2.15) and (3.16) are satisfied with c small enough and denote v := v — ul, then 
we get 



{u,DP)\\ep{T) < Cc\\dtUL, fiD'^UL, DPL\\^^^^,j,,^n/p-l^ 



+ C\\ul\\ , ,„ rr.nn/p. + C\\v\\ ^.on/ps \\UL 



/P\ 



L2{0,T;B;//) ' ^II''IIL2{0,T;B;7)I'-^HL2(0,T;B;/,^)- 
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This inequality together with the interpolation inequality 

irii II < ini"*"/^ irii"*"/^ 

I'''IIl2(0,T;B;//)II - ll^llLi(0,T;i?;//+^)"''"Loo(0,T;B;/^-l) 

ensures that $ maps {ul,\7 Pl) + Ep(T) (where Ep[T) is the "local" version of Ep) into itself 
whenever T satisfies 

Of course, for (3.16) to be satisfied, it suffices that to take R = c^/2 and to assume that T is 
so small as 

(3.20) ||Z)n,||^^(„^^^^„/.)<c/2. 

So if (3.19) and (3.20) are satisfied (conditions which depend only on the data) then one may 
conclude that ^ maps {ul,VPl) + E^{T) into itself. 

The proof of the contraction properties for $ in this context follows the same lines : we 
consider {ui,VPi) = '^{vi,VQi) with {vi,VQi) in {ul,VPl) + E^{T) for i = 1,2, then we 
bound all the terms Sfi, dgi and 6Ri as in the case of small initial velocity, except for 5Rs as it 
involves ||iii|L /„^ 6"/p-i\ which need not be small for T going to 0. For this latter term, we 

notice that, according to (A. 16), 

So we eventually get, 

\\{du,V6P)\\E,iT) < C(||l - Po||_^(^./P-i) + ||I)t;2||^^(o_^.^./P))||(<5u, V5P)|U^(T) 

+C'(/i-ii|(Si,VPi)|U^(T)+A^-^/=^||ni||^^(o^^^^n/p))||(*,V<X?)|U^(T)). 

Note that our assumptions on (ui , VPi ) and on the free solution ensures that if R and T 
have been chosen small enough then the factor of the last term is smaller than, say, 1/2. So 
the contraction mapping theorem applies. This completes the proof of the existence part of 
Theorem 2. Proving the stability and uniqueness follows from similar arguments. The details 
are left to the reader. 

3.5. Proof of Theorem 3. Given data {po,uq) satisfying the assumptions of Theorem 3, one 
may construct a global solution (it, VP) to System (1.2) in Ep. If Xu denotes the "flow" to u 
which is defined according to (1.9) then the results of the appendix ensure that, for all t € M+, 
Xu{t, •) is a C"*^ diffeomorphism of M". In particular, one may set 

p{t, •) := po o X-\t, •), P{t, •) := P{t, •) o X-\t, •), u{t, •) := u{t, •) o X'^t, •), 

and the algebraic relations that are derived in the appendix show that (p, u, VP) satisfies System 
(0.1). In addition, given that X„(t, •) is measure preserving and that DXu{t) — Id belongs to 

Bpf, the map a i-^ a o X^^{t) is continuous from B^^^ to itself if s G {n/p — l,n/p}, (see e.g. 
[8], Chap. 2). This implies that: 

• the Eulerian velocity u is in Cl)(M.-^-■, B"^ -^ ), 

• for any (/) G BpJ'^ , we have (j) ° ^t^it) G -^p,i ^ • So, given that po G M{B^^^^ ), we 
have 

Wn/p-l 



<pp{t) = {{<p o Xu{t))p^) o x-\t) G p;/i^ 

Hence p G L^{^+;M{B'''l'~^)). 
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• the chain rule ensures that 

VP=(^Au-VP)oX-\ 

So combining product laws and the invariance of B^i by right-composition, we get 

VPgLi(M+;^;//-'). 

• the chain rule also ensures that 

Vn= (^Au-Vu)oX-^. 

So using the fact that B"^i is a Banach algebra, and using invariance of B"^i by right- 
composition, we get Vti G Li(]R+; i?"j^^). 

In order to prove uniqueness, we consider two solutions {pi,ui, VPi) and {p2, U2, VP2) corre- 
sponding to the same data {po,uo), and perform the Lagrangian change of variable (pertaining 
to the flow of ui and U2 respectively. The obtained functions (ni,VPi) and (n2,VP2) both 
satisfy (1.2) with the same po and uq. Hence they coincide, as a consequence of the uniqueness 
part of Theorem 1. 

Appendix A. Appendix 

Let us first derive algebraic relations involving changes of coordinates. 

We are given a C^-diffeomorphism X over M". For H : R" -^ M™, we agree that H(y) = 
H{x) with X = X{y). With this convention, the chain rule writes 

(A.l) DyH{y) = D,H{X{y)) ■ DyX{y) with {DyX\, = dy^X\ 

or, denoting V^ = ^Dy, 

VyH{y) = {VyX{y))-V^H{X{y)). 
Hence we have 
(A.2) a,H{x) = DyH{y) ■ A{y) with A{y) = {DyX{y))-^ = D^X-^{x). 

Lemma 1. Let H be a vector-field over M". If we denote H = HoX then the following relation 
holds true: 

(A.3) diw^H{x) = d{Yy{A\DX\H){y) = dlw y{adi{DX)H){y) with x = X{y), 

where \DX\ stands for the determinant of DX, and adi{DX) for the adjugate of DX, that is 
the transpose of the cofactor matrix of DX. 

Proof: This stems from the following series of computations (based on integrations by parts 
and (A.2)) which hold for any scalar test function q: 

Jqix)di..Hix)dx = -lDM^).Hix)dx, 

D^q{X{y)) ■ H{X{y))\DyX{y)\dy 

Dyq{y)-{A\DyX\){y)-H{y)dy, 

q{y)diWy{A\DyX\H){y)dy. 

As ^ = \DyX\^^ a,(li{DyX) , we get the result. ■ 

Remark 1. Combining (A.2) and (A.3), we deduce that if a : M" -^ M then 



Axtt = divj^Vj^a = dw y{A\DX\\7 xa) = divy{a,di{DX) A\7ya). 
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Recall that if u is a time-dependent vector field with coefficients in Li(0, T; C^'^) then it has, 
by virtue of the Cauchy-Lipschitz theorem, a unique C^ flow Xy satisfying 

X,{t, y)=y+ [ v{t, X,{t, y)) dr for all t G [0, T), 
Jo 

and that X,^{t, •) is a C^'^-diffeomorphism over M". 

Lemma 1 enables us to deduce the following "magic" relation which is the corner stone of the 
proof of our main results: 

Corollary 2. Let v and w be two time- dependent vector fields with coefficients in Li(0,T;C^'^). 
Let Xy and X^ he the corresponding flows. Denote A^ := (DXy)^^ and A^ := {DX^^,)^^. Let 
us introduce the Lagrangian coordinates y^ and y^ pertaining to v and w, respectively, defined 
by 

X = Xy{y^) = Xy,{yu,). 
Assume in addition that 

\DXy\ = 1 and d\Y [A^w^) = with w^ := w o X^. 

Then \DXyj\ = 1 and for any C^ vector-field H, one has 

A\y H{x) = {DHy : Ay){yy) = div [A^H^^^y^) with H^ := H o Xy and H^ := H o X^,. 

Proof: With the above notation, the chain rule ensures that 

DxH{x) = Dy^Hy{yy) ■ Ay{yy). 

Hence taking the trace yields the left equality. 

Next, according to Lemma 1 and to our assumption over v and w, we have 

= div (Ayidy) = div(adj(Z)Xj,)tt;^) = divx'w{x). 

Hence Liouville theorem ensures that |DX^| = 1. So finally, applying Lemma 1 with X^ 
completes the proof. ■ 

Lemma 2. There exist n^ at least quadratic polynomials Pij : A4„(M) — > M o/ degree n — 1 
such that 

Id - adj(Id + C) = {C- (Tr C7)Id) + P2{C), 

where P2{C) is the n x n matrix with entries Pij{C). 

Proof: It suffices to use the fact that, by definition of the differential of adj, we have 

Id - adj (Id + C) = adj (Id) - adj (Id + C) = -dadj(Id)(C7) + P2{C). 

Now, 

adj (Id + C) = (Id + Cy^ det(Id + C) 

and the differential of the reciprocal operator at Id is C i— > —C while (idet(Id)(C) = (Tr C)Id. 
So d adj (Id) {C) = (Tr C) Id - C. ■ 

We now want to establish some a priori estimates for the flow which will be needed in our 
main results. The first difficulty that has to be faced is that when implementing the iterative 
process for solving (1.2), we are given the velocity field v in Lagrangian coordinates. Therefore, 
it first has to be checked whether the "flow" X„(t, •) deflned by 

(A.4) Xy{t,y) ■.= y+ v{T,y)dT 

Jo 

is a C^ diffeomorphism over M". This property is required for constructing the Eulerian vector- 
field V by setting v{t, ■) := v o X^^{t, •). 
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So let us assume that we are given some vector field v over [0, T) x M" with 

veai[0,T);B;;r'), dtveLii[0,T)-B;;r') and D^ € Li([0,T); ij;//). 
Differentiating (A. 4) with respect to the space variable yields 

(A.5) DX^it, y) := Id + / Dvir, y) dr. 

Jo 

As S"y'(M") is embedded in the set Co(^") of continuous functions going to at infinity, 
we deduce that X^ is a C^ function over M+ x M."". However, in general, Xy{t, •) need not be a 
C^-diffeomorphism over M" for all t € ]R_|-. So we assume that the smallness condition (2.5) is 
satisfied with c small enough. Then, using embedding we see that it guarantees that 

\\DX^{t, •) - Idlli^(Kn) < 1/2 for ah t G 1R+. 

Hence, for any t € M+, the map Xv{t, ■) is a local diffeomorphism. In order to show that it is a 
global diffeomorphism, we introduce the solution Yy to the ordinary differential equation 

(A.6) v{t,Y^{t,x)) + DX^{t,Y^{t,x))—Y^{t,x) = 0. 

at 

Under (2.5), the matrix DX^ is invertible at every point and (DX^)^^ — Id belongs to 

Looi^+'i^pi)- Indeed, one may write 

(DX^)-'^ - Id = ^(Id - DX^)''. 

fe>i 

Hence, using the assumptions over v and the product laws in Besov spaces (here we need that 
1 < p < 2n), (A.6) may be seen as an ordinary differential equation in the Banach space 
-^pi ' which may be solved on [0,T) according to Cauchy-Lipschitz theorem. In particular, 
differentiating Xy{t, •) o Yy(t, •) with respect to time, we easily gather that Xy{t, •) oYy{t, •) = Id 
for all t € [0, T). Therefore, one may eventually conclude that X^(t, •) is a C^ -diffeomorphism 
over M", with inverse l^(t, •). 

Let us now derive some "flow estimates" that will be needed for constructing the maps <5 and 
^. For completeness, the statements that we here state are slightly more general than needed : 
that \DX\ = 1 is not assumed. 



Lemma 3. Let p £ [l,-|-oo). Under Assumption (2.5) for v, we have 

(A.7) 
(A.8) 



(A.9) 

(A.IO) 

(A.ll) 



Id-adj(DX(t))|L„/p < \\Dv\\ ,^„/, 



at(adj(DX))(f)|L„/, < ||I)^(t)|L„/p, 

ai(adj(Z)X))(t)|| .„/,_! < \\Dv{t)\\B^/P-i if P< 2n, 
p,i p,i 

adiiDX{t)fA{t) -IdW^r./, <\\Dv\ 



^B^,P ^ ll^-llLi(0,i;B;//)- 
Proof: According to Lemma 2 and to (A.5), one may write 

ld-ad]{DX{t))= I {Dv-diYvld)dT + P2(( I Dv dr 

where the coefficients of P2 are at least quadratic polynomials of degree n — 1. Given that B'^-^ 
is a Banach algebra and that (2.5) holds, we readily get the result. 
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In order to prove the second estimate, we just use the fact that, under assumption (2.5), we 
have 

(A.12) A{t) = (Id + Cit))-^ = ^(-l)'=(C(t))'= with C(t) = / DvdT, 



keN 



and that B^i is a Banach algebra. 

In order to prove the third inequahty, we use the fact that, according to Lemma 2, we have 

dt (adj (DX)) = — ( Id + / (div t; Id - Dv) dr + Pil f Dv dr 



Hence 



dt{a.d]{DX)){t) = {d\wv{t)ld-Dv{t)) + dP2( Dvdr) ■ Dv{t). 
As the coefficients of dP2 are polynomials of n? variables that vanish at 0, we get 



\dtiadiiDX))it)\\^,, < \\Dv{t)\\^„/,(l + 



t 
Dvdr 








hence (A. 9). Proving (A. 10) is similar: it is only a matter of using the continuity of the product 
from B;lf-' X ij;(^ to B;lf-\ if p < 2n. 

For proving the last inequality, we use the decomposition 

adi{DXfA - Id = {adi{DX) - IdfA + '^{A - Id). 

So combining Inequalities (A. 7) and (A. 8), and the fact that B"^i is a Banach algebra, we get 
the result. ■ 

Lemma 4. Let vi and V2 he two vector-fields satisfying (2.5), and 6v := 'U2 — t^i. Then we have 
for all p G [1, +oo), 

(A.14) \\adi{DX2) - adj(I)X,)||,^(^^^^./.) < \\DH,^^^^.^^;f.y 

(A.15) \\dtiadi{DX2) - adi{DX,mL,(^^^.^B;{n - \\^^kr{R+;B;{y 

(A.16) \\dt{adi{DX2) - adj(Z)Xi))||^^^j^^.^„/,-i^ < ||D&||^^^^^.^„/,-i^ if p < 2n. 

Proof: In order to prove the first inequality, we use the fact that, for i = 1, 2, we have 

Ai = (Id + d)-^ = V(-l)^Cf with Ci{t) = [ Dvi dr. 
k>o ■'^ 

Hence 

a2-a^=y: (ci - ci) = (fD5v dr) y: e CiCf-'. 

k>l ^''^ ^ k>lj=0 

So using the fact that B^^f is a Banach algebra, it is easy to conclude to (A. 13). 
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The second inequality is a consequence of Lemma 2 and of the Taylor formula which ensures 
that, denoting dC := C2 — Ci, 

adj(Z)X2) - adj(Z)Xi) = (Tr(«7))Id - dC + dP2{Ci){dO) + ]^(fP2{Ci,Ci){dO,dO) + • • • 

where the coefficients of P2 are polynomials of degree n — 1. As the sum is finite and B^i is a 
Banach algebra, we get (A. 14). 

In order to prove the last two estimates, it is only a matter of differentiating the above relation 
with respect to t. Keeping in mind the definition of Ci and C2, we get 

9i(adj(L>X2) - adj(Z)Xi)) = (div(fc)Id - DSu + dP2{Ci) ■ dtdC + (fP2{Ci){dtCi,dC) + ■■■ . 

Then using the product laws in Besov spaces yields the desired inequalities. ■ 

Finally, we have to justify that the multiplier space M.[B^'^ ) contains characteristic func- 
tions of C^ bounded domains, if p > n — 1. This is a consequence of the following lemma. 

Lemma 5. Let Q. he the half-space M" or a bounded domain of M" with C^ boundary. Assume 
that s € M and p,q (z [1, 00] are such that 

(A.17) -l + -<s<-- 

p p 

Then the characteristic function xn of Q belongs to the space A4{Bp JW"")). 

Proof: This result which belongs to the mathematical folklore is closely related to the fact 
that under Condition (A.17), functions in Bp (Q) extended by on the whole space, belong to 
Bp (W^). In the case where Q is the half-space M" the lemma has been proved in [7], Prop. 3. 
If fi is a bounded C^ domain, then one may find a finite number N of C^(M") functions (pi 
and C^ diffeomorphisms ipi so that for any u £ Bp JW^), 

N 

uln = ^ ucpiln and (ucpiln) o ipi = l^+ ■ {{u<pi) o ip^) . 

i=\ 

Now, because Condition (A.17) is satisfied, the space i?p„(IR"') is stable by multiplication by 
smooth compactly supported functions, and by C^ change of variables (see e.g. [8], Chap. 2). 
Therefore the functions {u^i) o t/^j belong to 5^ (M"), too. Using again the stability of this 
space by multiplication by 1ir+, one may thus conclude that (u^jlf^) o ip^ £ Bp^g{M.'^). Hence 
upiln is in B^ „(IR") for all i G {1, • • • , N}. This completes the proof of the lemma. ■ 
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